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Abstract 

We propose to minimise the cosmic confusion between Gaussian and non Gaussian 
theories by investigating the structure in the m's for each multipole of the cosmic 
radiation temperature anisotropics. We prove that Gaussian theories are (nearly) the 
only theories which treat all the m's equally. Hence we introduce a set of invariant 
measures of "m-preference" to be seen as non-Gaussianity indicators. We then derive 
the distribution function for the quadrupole "m-preference" measure in Gaussian 
theories. A class of physically motivated toy non Gaussian theories is introduced 
as an example. We show how the quadrupole m-structure is crucial in reducing the 
confusion between these theories and Gaussian theories. 



Revised Version 



1 Introduction 



Several competing theories of the Early Universe have now been proposed. Overall, two 
major families have emerged, one centred around the idea of cosmic inflation [Q], the other 
relying on topological defects formed in phase transitions 0]. Crudely speaking, it was 
found that the first gives rise to a Universe with Gaussian features, whereas the latter 
should leave some sort of non-Gaussian imprint in the Universe. Producing a cosmic 
crucible has however been problematic. The discovery of temperature fluctuations in the 
cosmic microwave background (CMBR) opened doors to a new field with a chance for a 
decisive experiment. In this letter we examine the obstacles cosmic variance raises to the 
design of such an experiment and we propose a new approach to circumventing 
them. 

Cosmic variance in the CMBR comes about because for the Universes predicted by 
every theory (and probably for our Universe too) there is only a probability distribution 
function for what the CMBR sky should look like. Different observers, located at different, 
uncorrelated, points of the Universe, observe different skies. These follow a distribution 
function which varies from theory to theory, but which is never so peaked that only one 
sky is possible. The distributions naturally overlap, and so cosmic confusion arises. Mere 
mortals, with access to only one sky, may be left with the problem that their sky could have 
derived from totally different theories. To be quantitative, let us define cosmic confusion as 
the percentage of common skies generated by two theories. Then, if a measurable quantity 
Q is predicted to have distributions Fi{Q) and F2{Q) in the two theories Ti and T2, the 
cosmic confusion in Q between Ti and T2 is 

Cq(Ti,T2) = / dQ min(Fi(g),F2(g)) . (1) 

Jail Q 

C(3(Ti, T2) varies between (no doubts) and 1 (totally confused), and is the confusion before 
Q is actually measured. It is the quantity to be minimised (by means of an appropriate 
choice of Q) when projecting experiments which, one hopes, will be conclusive. Once Q 
is measured, say, with an outcome Q G Vq, Vq = {Qq — AQ~ , Qq + AQ"*"), the confusion 
becomes 

and if /-p dQ Fi > dQ F2, the probability of Ti over T2 is 

Iv dQ Fi 

The attitude adopted in this letter is alarmingly anti- anarchist. We seek to minimise the 
cosmic confusion between Gaussian and non Gaussian theories using quantities associated 
with the low order multipoles of ^ (as opposed to / > 30, where confusion is known to be 
small again). Rather than looking at the angular power spectrum or the quadrupole 
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intensity Qrms, we investigate the structure in the m's for each / (we will actually concen- 
trate on / = 2 in this letter). Gaussian theories are known to treat all m's in the same way. 
In Section |^ we prove that such a feature is in fact peculiar to Gaussian theories, within a 
large class of theories. The physically most relevant theories, however, appear to be outside 
the class considered, but in Section ^ we show how "m-preference" still plays an impor- 
tant role as non-Gaussianity indicator in the general case. Hence quantities measuring 
m-preference appear as good candidates in our quest for low cosmic confusion. In Section 
^ we write down the most general multipole invariants and identify a set of invariant mea- 
sures of m-preference. We then speciahse to the quadrupole (Section H) and in Section ^ 
the distribution function for its measure of m-structure is derived for Gaussian theories. 
A class of physically motivated toy non Gaussian theories is introduced in Section ^ We 
exhibit one extreme case where the confusion in Qrms is !> but where the confusion in the 
quadrupole m-structure is 0. We conclude with a brief summary of the results obtained. 



2 Structure in the m's as a sign of non Gaussianity 

In analyzing ^(6*, 0) maps it is traditional to use the multipole expansion: 

STT-I OO ( OO / 

^ = E E «!n>l = E E bLQL. (4) 

1=0 m=—l 1=0 m=—l 

The real and complex spherical harmonics (Q^ ^'^^ ^m) related by 

yo' = Ql ; yL = (T)'" ^^^^-- (5) 

and similar relations apply to their coefficients and a^. In this Section we assume 
that the 6^ are independent random variables. Under this condition we prove a rather 
surprising result: that only Gaussian theories are statistically spherically symmetric (SSS) 
and m-amorphous (i.e: treat all the in the same way for each /). The implication is 
that conversely any non-Gaussian theory is to some extent anisotropic, favouring particular 
directions in the sky and some m's over the others. 

We say that a theory is SSS if any two quantities in the sky related by an 0(3) trans- 
formation are predicted to have the same distribution function {f{R{Q)) = f{Q)). We 
also say that a theory is m-amorphous if for a given fixed / the 6j„ are independent random 
variables with the same distribution function, regardless of m and the choice of axes. If 
this distribution is Gaussian we then say that the theory is Gaussian. Let us first prove 
the following Lemma: 

Lemma 1 Let x and y be two independent continuous random variables with the same 
distribution f . Let z = ax + j3y with + (3"^ = 1. Then, if f is also the z distribution, f 
must be a Gaussian. 
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g{t)=Y.grr. (8) 
(?„ = (a" + /3"K (9) 



Proof : Let be the characteristic function of / (e.g. [Q). Then 

= <P{at)<P{f3t) . (6) 

Putting g = log / this imphes that 

g{t) = g{at) + g{(3t) . (7) 

Now expand g in power series: 
Then (|^ requires that 
and so (7„ 7^ iff 

a" + /3" = l. (10) 

But we know that + = 1, and by studying the curves ( pUj ) in the (a, P) plane for all 
the n's one finds that they only intersect at a = or /3 = 0. Hence n = 2, and so g{t) oc 
and oc e^^* , which is the characteristic of a Gaussian distribution. (Q.E.D.) 
We can now prove the central result of this Section. 

Theorem 1 For a theory in which all the 6^ are continuous independent random variables 
the following three properties are equivalent: SSS, m-amorphism, and Gaussianity. 

Proof: It suffices to show that SSS implies Gaussianity. (It looks as if we are proving 
more than the equivalence of m-amorphism and Gaussianity because it is not immediately 
obvious that a SSS theory is also m-amorphous. This is because given a basis {Ql^} 
its elements cannot all be 0(3)-transformed into each other.) We start by noting that by 
rotating an angle around z, for a given fixed m 7^ 0, the elements and transform 
like 

Qm ^ Qm ^Qm PQ—in 

Qim - QL^ = f3Ql + aQL^ (11) 

with a = COS0, and (3 = sin (p. Choosing = 7r/2 we then know that f{bl^) = /(&'_„) = /. 
Since also f{bl^) = f for any 0, using the Lemma we know that / has to be a Gaussian. 
Now, let us apply the most general rotation to Qj^. Any element of its orbit will have the 
distribution function /. The subspace spanned by the orbit is not the whole space, but 
because the representation is irreducible, it will not be an invariant subspace. This means 
that vectors V outside the subspace exist which, although not obtainable directly by a 
rotation of Q''^, can still be obtained by rotating a linear combination of rotations of Q^- 

V = R[Y^aMQl)] . (12) 

n 

By choosing = 1 one finds that fiV) = f. V and the orbit of Ql^ still do not span 
the whole space, but the procedure described has allowed us to extend non trivially the 
subspace spanned by the elements for which the distribution is /. By iteration one can 
therefore show that for all m the distribution /(&^) is a Gaussian distribution function / 
with a variance which can only depend on /.(Q.E.D.) 
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3 SSS non-Gaussian theories 



The theories for which all the 6^ are independent form a large class. Within this class m- 
structure is a certain indicator of non-Gaussianity. However only non-SSS non-Gaussian 
theories can belong to this class. While violations of SSS have appeared in proposed 
models (such as the ^ brought about by any anisotropic cosmological model), one should 
beware of these violations as they can never be tested experimentally. Furthermore most 
practical applications concern ^ due to perturbations in FRW models, for which SSS is 
always satisfied. If one takes SSS as a starting point, then one must confront the following 
corollary of Theorem 1. 

Corollary 1 For SSS non- Gaussian theories the 6^ can never be all independent random 
variables. 

Once one allows the 6^ to be dependent the rest of Theorem 1 breaks down (see the end 
of Section |^ for an example). For SSS theories m-structure may or may not be associated 
with non-Gaussianity. In some cases m-structure is still the hallmark of non-Gaussianity 
(see Section |^), but there is also a class of non-Gaussian m-amorphous theories. 

However we now show that if a non-Gaussian theory is m-amorphous and sufficiently 
distant in functional space from Gaussian theories, then the C' must be distinctly non- 
Gaussian. The power spectrum appears to save the day when m-structure is useless. We 
show this fact with the aid of the concept of usefulness of a variable. Let the distance 
between two theories be 

V{T,,T,)= I IF1-F2I e(0,2). (13) 

J all variables 

For any subset of variables S we have 

i-p/2<C5(ri,r2)<i, (14) 

showing that the minimal possible confusion is Cmin = 1 — ^^/2. If P = 0, the two theories 
are the same, and no wonder we are totally confused whatever set of variables we look at. 
For X> 7^ 0, however, one can assess the usefulness of the variables S with the quantity 

7 / 1 ~ Cmin i^r-\ 

Us = l- ——, . (15) 

Us varies between (utterly useless) to 1 (perfect choice). Obviously if S contains all the 
variables. Us = 1, but what we should be after is maximising U with the smallest number 
of variables. Now if Ti and T2 are two m-amorphous theories (Gaussian or not) we have 
Fibi) = fiC^) and so 

1^ = J (UdbLmibl) - F,ibl)\ = J {Y[dC')\F,{&) - F,{C')\ . (16) 

Im I 

Therefore Uqi = 1, implying that if Ti is Gaussian and ©(Ti, T2) ~ 2 then Cc!(Ti, T2) ~ 0, 
that is, C' for T2 must be distinctly non-Gaussian. 
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The important converse statement of the last paragraph is that a very non-Gaussian 
theory with a very Gaussian must display a distinctly non-Gaussian m-structure. This 
clarifies the role of m-structure for SSS theories: m-structure may or may not work for 
them as a non-Gaussianity indicator, but it always works when does not. 



4 Multipole invariants and structure in the m's 

These results prompt us to quantify m-preference as m-preference measures will also mea- 
sure non-Gaussianity. For this purpose it will be useful to first recall the isomorphism 
between the space spanned by the real spherical harmonics of degree / and the space of 
the real traceless symmetric cartesian tensors of rank / {Qi...j = Q{i...j), Qiij...k = 0). This 
isomorphism can best be established by writing the functions Ql^{6,(f)) in cartesian coor- 
dinates. The resulting polynomial of degree / has the form Q["^jX^...x^ , identifying a set of 
traceless symmetric tensors of rank /: 

(17) 

The isomorphism can then be extended to the rest of the space using linearity. The basis 
{Q^i^j} is known as an irreducible tensor set and one can check that indeed both spaces 
are 21 + 1 dimensional (see for instance 0). 

Now for whatever combination of the (or bl^) one considers it is important to demand 
spherical invariance, so as to factor out the artifacts introduced by our inevitable choice 
of axes. Hence we look for multipole invariants. Under an axes transformation defined by 
the Euler angles [ijj, 6, (p) the a^^ transform according to 

"'m / y m'm m' 



where D^^i^^ip^ 6, (j)) is the Wigner matrix. The tensors Qi...j, on the other hand, transform 
just like any other cartesian tensor. The fact that the matrices -D^'m form an irreducible 
representation of the 3-dimensional rotation group requires, by means of Schur's Lemma, 
that the only matrices they commute with be multiples of the identity. As a result, the 
only independent bilinear invariant one can construct out of the a^^ (or tf^) is 

c^= i: |«Lf = E (bir (19) 

m=—l m=—l 

more proverbially known as the angular power spectrum. However, by simply counting 
degrees of freedom, one is easily convinced that there are in fact 2/ — 2 invariants for each 
I >2. What Schur's Lemma states is that 2/ — 3 of them will not be bilinear (that is, they 
will be of higher order). Therefore, the invariants are in general multilinear forms 

21+1 
i...j=l 
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^i'.lj' — ■ (21) 



satisfying the condition 

j(n 

i' ...j' I'f f j-^i...] 

These are very difficult to find for a general /. Since the rotation matrices -D^'m ^^^^ ^^^o 
a subgroup of SO {21 + 1) we know that / of these invariants will be the fundamental 
representations of the Casimir operators of S0{21 + 1). Still, even these have a rather 
unpleasant form. We have found it easier to write the invariants in terms of the tensors 
{Qi...j} instead. They then appear as the independent contractions of a symmetric traceless 
tensor of rank / in 3 dimensions. We cannot produce a general formula for these, but this 
approach has proved to be more manageable case by case. 

What is the meaning of the The invariant is obviously a measure of the overall 
intensity of the multipole. It treats all m's equally, so it does not assess m-preference at all. 
The other J^"-*, however, treat the various m's differently (see Section ^ for an example). 
For all invariants of order m > 2 one should then define the ratios: 

(n) 

The shape factors r^"'' reveal how some m's are preferred over others, or, in other words, 
how much and what type of anisotropy there is in the multipole, but only in so far as there 
is a rotationally invariant meaning to the concept. Demanding rotational invariance is 
important since otherwise our anisotropy measures would refiect not only the directionality 
pertaining to the CMBR but also the directionality imparted by our choice of axes. To 
make the point clear consider the case / = 1. Using (|1^) the dipole can be seen as a vector, 
for which the only invariant is (the dipole has no shape). Clearly, any dipole has a 
direction, so it could never distinguish between an isotropic and an anisotropic theory. 

There are also invariants which combine different /'s. By counting degrees of freedom 
we find that, for each pair /, I' > 2, there are 3 inter-Z invariants which do not depend on 
the invariants for each /. We may wish to look at J*^"^ as generalised eigenvalues. For a 
general I one can define the eigenvectors of Qi,„j as the set of axes where three independent 
prechosen components of Qi,,,j are set to zero (for I = 2, Qij = ior i ^ j), and its 
eigenvalues as the values of the remaining components in that basis (for / = 2 the diagonal 
components). The inter-Z invariants can then be interpreted as the Euler angles of one 
set of eigenvectors with respect to the other. They are uniformly distributed in Gaussian 
theories, but not in defect theories, where they reveal the correlations between successive 
generations of defects P]. 



5 Quadrupole variables and invariants 

We now specialise to / = 2. Using the isomorphism (|r^) a quadrupole can be seen as a 
traceless symmetric real matrix Qij. This can be fully characterised by the Euler angles 
{ip,9,(f)) specifying the rotation required to diagonalize Qij, together with its diagonal 
form Qij. The matrix Qij has two invariants under 0(3) transformations, which can be 
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parameterized in a variety of ways. The eigenvalues Qij = diag(Ai, A2, A3) (subject to 
Ai + A2 + A3 = 0), if considered in modulus and in unordered triplets, constitute one such 
parameterization (recall that improper rotations permutate the eigenvalues). One can also 
expand Qij in multipoles: 

Q,, = aQl + bQl. (23) 
Again, a and b are not left unchanged by improper rotations, but transform according to 



(1) 
(2) 

(3) 



a 
b 



—a 
-b 
VSb 



2 

■\/3a + b 



(24) 



Modulo these transformations, a and b also parameterize the Qij invariants. Finally, one 
can do away with any residual transformations by considering the invariants: 



h = TtQ' 



lTrg3 



: IdetQI =nA. 



U67r 



30 



[a^ + b^ 



2bo 

61 -f,+b. 



(25) 



b-2 



b-1 
b-2 

bo h 



/ 15 ^ 


J 3/2 




U67r. 




Vs\3 J 



(26) 



We should note that since in a SSS theory {Aj} and (a, 6) are subjected to residual sym- 
metries interconnecting them they can never be independent random variables. In fact 
their joint distribution functions, subjected to the same symmetries, can never factorize 
(cf. (^3])). Therefore, for statistical applications, they are not good parameterizations of 
the quadrupole invariants. Also, one can prove that 



A1A2A3 



(A? + Al + Ai) 



2\3/2 



< 



3^6 



and so we have the constraint 



1 2 < 



3/2 



3^6 ■ 



(27) 



f28) 



Consequently, even though Ji and I2 are not interconnected by residual symmetries (and 

so their joint distribution function does factorize), their ranges of variation are dependent 

j-3/2 

(Ji e (0,00), I2 e (0, 1^)). Hence h and h are in fact dependent random variables, a 
fact manifest in that the integrated distribution function of I2, F{l2) = J F{Ii, l2)dli, is 
different from the factor ^2(12) appearing in F(/i,/2) = Fi{Ii)F2{l2)- For this reason we 
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shall characterize the quadrupole Qij by the Euler angles (■?/', 6*, 0), the intensity Ji (known 
as Qrms)) and the ratio r (to be called the "quadrupole shape"): 

r = 3^6-1^ (29) 
-'1 

which varies in the range r G (0, 1). These five quantities are statistically independent: 
their joint distribution function factorizes and their ranges of variation are independent. 
The shape factor r is an invariant measure of how axis-symmetric the quadrupole is. 
It varies from maximal symmetry for r = 1 (two equal A; 6 = possible) to maximal 
symmetry breaking for r = (a null A; a = possible). 



6 The distribution of the quadrupole shape in Gaus- 
sian theories 

It turns out to be easier to find first the distribution function -F(a, b,tlj,6,(j)). Since Gaussian 
theories are SSS we know that F has to have the form 



F{a,b,ip,9,(l)) = F{a,b) 



sin 6 



(30) 



where F{a, b) is further invariant under the transformations (p4D. Then start by writing: 

d{bo,bi,b^i,b2,b_2) 



F{a,b,^,e,(j)) 



F(6o,&l,&-l,&2,&-2 



(27r)5/Vi ■ 



d{a,b,^/j,e,(l)) 
d(ao, ai, a_i, 02, a_2 



d{a,b,'^,e,(j)) 



(31) 



Now let D'^^,{ilj,6,(f)) be the Wigner matrix of the rotation taking Qij to Qij. Then 



= aD'^Q + bH^ where = (-0^2 =t ^In--!)!^- The matrices -Dmm' can be seen 



as generalised spherical harmonics satisfying differential equations similar to the spherical 
harmonics equations In particular: 



dil) 
86 



mcot eDi^Q + J{2-m+l){2 + m).e"^.D. 



m-l 



(32) 



and 



-iniH. 



2± 
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m 

dip 

m 

de 



2cos dHi+ - 



^ 0-2- 



These relations imply that the determinant in (^Tj) must have the form 
d{ao, ai, a_i, 02, a-2) 



d{a,b,ip,e,(f)) 



(33) 



(34) 



and its exact form can then be found from the symmetries. In fact 
01) 02, 03 oc sin6' and (0) leads to: 



F(a,6,V,^,0) 



„ sin 61 



a^ + b^' 



b{3a^ - b^ 



]) implies that 



(35) 



Having found F{a, b, ip, 6, 0) one can now write F{Ii, I2, ip, 0) as 



F(Ji,/2,^,^,0) 



^ F(a,6,V^,^,0) 

8 branches 



(9(a, b) 



(36) 



The branches referred to under the summation sign are the 8 branches generated by the 
symmetries (|2^ which give the same values (Ji,/2). It can be proved that not only 
F{a,b,ip,6,(p) but also the determinant 



d{a,b) 



are invariant under (|2 



r \ xiivuxxuiii. ^^^yj.^^ Hence all 8 

d(Ii,h) ^—^^ 

branches give the same result which can be more easily evaluated for the a > \^b > 
branch. One then gets 

F(/i,/2,V^,5,0) = C'^.e"4 (37) 



871' 



with (72 = (T2\/-^, and from it one finally obtains 



Wit 



2ct, 



e '*2Ji/^ sin^ 



F{h,r,iP,e,< 



(38) 



where the proportionality constant was computed from a normalization condition. We 
have thus determined the distribution of the 5 independent variables of the quadrupole 
in a Gaussian theory. The quadrupole intensity Ji is a xl variable (as well known), its 
shape r is uniformly distributed in (0, 1) and its axes variables {ip, 6, 0) have a uniform 
distribution. We have confirmed these results with a Monte-Carlo simulation. 

We can now, as a side remark, provide the example required in Section |^. Note that for 
/ = 2, SSS (equivalent to F(/i, r, ^/^, ^, 0) = /(/i,r)|s|) implies that F(6o, &2, ^-2) 

3 /2 

oc f{Ii,r)/Ii . On the other hand, m-amorphism (equivalent to F(6o, &2, ^-2) = 

/(/i)) implies that F{Ii,r,ip,6,(p) oc /(/i)/f^^ sin These general formulae show that, 
unless we require that F(6o, &i, ^2, ^-^2) factorizes, SSS does not imply m-amorphism, 
and neither implies Gaussianity. 
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7 A (toy) non-Gaussian theory 



We now produce a toy non-Gaussian theory in order to show how m-structure measures can 
minimize the confusion between Gaussian and non-Gaussian theories. Imagine a theory 
which imprints a quadrupole in ^ such that there is always an eigenvector basis where 

Qij — ciQij 5 (39) 

but the eigenvectors' orientation is uniformly distributed. This toy model can be physically 
motivated. It is a very good approximation to the quadrupole brought about by a B-field 
with a very large coherence length (see or fill)- It is also a good approximation to the 
quadrupole predicted in a texture scenario with scaling and with a small average number 



of defects per horizon volume {< N > <^ 1). Reasoning in accordance with |12] one should 
expect the quadrupole in this scenario to result mostly from the gravitational effects of the 
last defect. If most of the collapses are reasonably spherical the quadrupole will then be 
approximated by (PU]). 

The quadrupole invariants for this theory are 

h = —a' 
IGtt 

r = 1 . (40) 

Consider then the extreme situation in which a is such that Ii is x|-distributed. The other 
quadrupole variables follow the distributions 

/(r) = 25(r-l) 

me,<P) = f^. (41) 

Hence, the cosmic confusion in the quadrupole invariants between this theory (T2) and a 
Gaussian theory (Ti) with the same Ji is 

C/,(Ti,T2) = 1 

a(Ti,T2) = 0. (42) 

The distributions fTii^n) are somewhat complicated (they will be derived and used in p|), 
but it should be obvious that they do not reduce the confusion to 0. The invariant r is 
clearly the appropriate variable to thoroughly make out the difference between the two 
theories. Using the concept of usefulness Uci = and Ur = I. 

Naturally one should not expect to find such an extreme situation, say, in the texture 
scenario. First, the confusion in Ii is probably not as big, and second r will peak around 
1 but not so drastically. The current all-sky simulations have not yet produced a sufficient 
number of skies for distribution functions to be predicted with any statistical relevance 
(only 8 skies were produced in 0). In p we shall describe a computational shortcut 
which will allow a prediction for Ftex{Ii) and Ftex{f^)- 
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8 Summary and planned developments 



In this letter we proposed yet another non-Gaussianity indicator: m-structure (see [jT3 
for a summary of other options). We gave general arguments linking m-structure and 
non-Gaussianity, and set up the general framework for quantifying m-structure. We wrote 
down explicit measures for the quadrupole m-structure, and derived their distribution in 
Gaussian theories. We finally displayed physically motivated examples of theories for which 
only m-structure can bring out non-Gaussian nature. 

The signal to noise ratio in current all-sky data is not yet good enough to allow any 
meaningful measurement of r (still we hope to apply our approach to COBE data analysis). 
Nevertheless, an accurate measurement of r meets no obstacles of principle, with the 
possible exclusion of the effects of galactic obscuration. It is not clear how much of a 
problem this is, but we hope to generalise our work to /> 2, so as to minimise this effect. 

In work currently in progress we also seek to derive an approximation to /(r) in var- 
ious defect theories. Preliminary results seem to generally link the quadrupole with one 
individual defect (the closest defect). That being the case, the defect morphology has an 
imprint on the quadrupole shape, as the defect symmetries impose selection rules on the 
spherical harmonics expansion. A rich peak structure in /(r), depending on the defect 
type, ensues. This supports the claim that r, if measurable, is a low confusion variable 
when confronting defect theories among themselves and against Gaussian theories. 
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